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Abstract 

We present a theoretical method for the design and optimization of quantum corrals with spe- 
cific electronic properties. Taking advantage that spins are subject to a RKKY interaction that 
is directly controlled by the scattering of the quantum corral, we design corral structures that 
reproduce spin Hamiltonians with coupling constants determined a priori. We solve exactly the 
two-dimensional electron gas scattering problem for each corral configuration within the effective 
mass approximation and s-wave scattering using a Green function method. Subsequently, the ge- 
ometry of the quantum corral is optimized with an algorithm that combines simulated annealing 
and genetic approaches. We demonstrate that it is possible to automatically design quantum cor- 
rals with complicated target electronic properties, such as multiple mirages with predefined relative 
intensities at specific locations. In addition we design structures that are particularly sensitive to 
the phase shift of impurities at certain positions allowing the measurement of the value of this 
parameter on the copper surface. 
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I. INTRODUCTION 



The fabrication of solid state quantum computing devices involves both the ability to ma- 
nipulate mater at the nanoscale l| and to design structures that interact according to pre- 
scribed quantum computing Hamiltonians^j. Tuning and controlling quantum interactions 
are basic problems faced by the quantum computing community. This requires engineering 
the structure of electronic wave functions in different environments and conditions. The 
manipulation of individual atoms with STM has made possible the construction of arbitrary 
quantum structures on top of surfaces. In particular, after the seminal work of Crommie 
et al.ljsl, it has been possible to build quantum corrals of arbitrary shape, placing atoms 
one at a time. Quantum corrals are a collection of atoms arranged in a controlled manner 
on top of a metallic surface. These novel structures generate quantum confinement of the 
surface conduction electron wave functions leading to striking phenomena such as resonant 
electronic states and the formation of quantum mirages [4J. Generally speaking, quantum 
mirages are the projection of a perturbation on a point into another distant point of the 
surface. Manoharan et al. ^ were able to partially project the Kondo P| electronic structure 
up to II2A away from the actual Kondo impurity. In these experiments, the formation of a 
Kondo mirage relies on the focusing properties of elliptical corrals. These findings triggered 
a collection of theoretical papers modeling and extending the concept of quantum 

,„.agesBQ. 

n 

This new phenomenon has been proposed as a tool for remote probing 14^ and as a 

n □ 

way to enhance j8|, |9|| the interaction between localized magnetic impurities on a surface. 
Unfortunately, the design of quantum corrals with certain desired properties is still a matter 
of trial and error j^. We think that both technological applications and novel physical 
phenomena could emerge if one could tailor the electronic wave functions that are responsible 
for the quantum mirage formation. For that purpose we believe that new systematic methods 
to design quantum corrals are required. 

In previous works it was shown that the magnetic interaction between magnetic impurities 
can be strongly enhanced due to the electronic confinement produced by quantum corrals ^, 
0]. Our present objective is to generate the optimal geometrical arrangement of quantum 
corral atoms that achieves an arbitrary predefined quantum mechanical interaction. This 
type of problem is usually denoted as inverse problem. In contrast to the direct problem, 
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where a physical property is the result of the geometry of the material, in the inverse problem 
the optimized geometry is a result of the physical property targeted. Inverse problems have 



been addressed in a variety of systems ranging from band gaps of solids design of 
antennas 12], or intermolecular potential fitting 13]. In order to solve the inverse problem one 
needs two ingredients: i) a fast method to solve the direct problem (i.e. electron gas response 
for a given corral configuration) ii) an efficient algorithm to select new configurations and 
optimize the structure to achieve the target properties. 

In this paper we address the 'inverse quantum corral problem'. We concentrate in the 
search of structures that could reproduce predefined magnetic Hamiltonians with desired 
coupling constants by means of an enhanced RKKY interaction. We demonstrate the pos- 
sibility to physically construct several Hamiltonian examples. 

Our theoretical approach consists of i) solving each trial corral configuration with a Green 
function method]^ and ii) optimi zing the corral geometry with a combination of simulated 
annealing and genetic algorithms We find that i) it is possible to design quantum 
corrals with multiple mirages with both predefined positions and intensities once one has an 
accurate model of the scattering, ii) Our theory allow us to design quantum corrals that are 
particularly sensitive to the values of the phase shift, the only free parameter of the model. 
We show that it is possible to design structures that can be used to measure this phase shift. 

II. DIRECT PROBLEM SOLUTION 
A. Solving the scattering 

We consider quantum corrals on the Cu(lll) surface. Cu atoms on this surface form a 
triangular lattice with spacing a = 2.55A. This underlying triangular lattice defines also 
a triangular lattice for the equilibrium positions of the impurities that form the corral. 
The electronic structure of Cu(lll) consist of a surface band with a Fermi energy ep laying 
450 meV above the bottom of this band. We will treat the surface band as a two-dimensional- 
electron gas in the effective mass approximation. There is broad consensus 
llfil that the quantum mirages form as a result of the interplay between the impurities at 
the surface and this two dimensional surface band. 

Several types of impurity atoms have been used to construct quantum corrals, examples 
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are Fe P| and Co 0| . In all cases, the Fermi wave length of the electrons in the surface band 

is much larger than the spatial extension of the perturbation introduced by a single impurity 

atom. Therefore, an s-wave scattering approximation is good enough for the description of 

I I 

the electronic scattering by the impurities 16]. Accordingly, the impurities that form the 
corral are described by a single parameter being the s-wave phase shift 6. 

For the case of scatters located at positions |a,;| the retarded Green function C'^* of 



the system can be expressed in the following way 



17|: 



G-Hr,r') = Gf(r,rO+ (1) 
+ X^Gf(r,a,)[lA-G;,]r;Gf(a„r') 

where [Ggjij = (1 — 6ij) x G'Q'^*(aj, a^) and t is the t-matrix of the identical scatters; in 
terms of the phase shift 6, t = i^(e*^^ — 1), being m* the effective mass (m* = 0.38me). 
In two dimensions and within the effective mass approximation G}^^{r,r';u) = 

^(Yo{k\r — r'\) — 2Jo(fc|r — r'|)), hk = \/2m*uj. While this method can be equally effi- 
cient for any phase shift model, for the purpose of this study, we chose to set 6 = ioo (unless 
otherwise specified) which has been determined to produce the best fit of the experimental 
information for quantum corrals constructed with Fe atoms on the cooper surface jl^. 

Equation (P) provides a straightforward way to obtain the exact Green function after 
an X matrix inversion. For the typical number of scatters that can be experimen- 
tally handled (A^ ~ 90) this method is fast enough to solve thousand of different atomic 
configurations in seconds, each of them representing one trial step. 

The Green function contains all the information about the solution of the independent 
particle problem and can be used to compute, for example, the local density of states 
(LDOS): 

V{r;uj) = --lmG''\r,r-ij) (2) 

or the perturbation caused by additional atoms, in particular magnetic impurities. 

Let us note that G^'^^ is a function of the atomic positions {aj}. In this work we address 
the question of which is the corral (i.e. the set {aj}) that gives predefined couplings between 
magnetic impurities by means of modulating the RKKY interaction. 
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B. RKKY interaction in two dimensions 



The exchange interaction of a magnetic impurity at coordinate Ri with an electron gas 
with couphng constant Ji is described by the Hamiltonian: 

i^imp = -JiSi • V'^(Ri)o-V'(Ri) (3) 

where Si describes the impurity spin, cr = (cr^;, cr^, cr^) are the Pauh matrices and V'^(r) = 
('0|(r), '0j(r)) describes the electron field. 

If a second magnetic impurity with spin S2 is placed at R2, the effective impurity- impurity 
coupling mediated by the electron gas can be written as: 

H12 = - J12S1 ■ S2 (4) 

with 

J12 j du;/(a;)Im[G'f (Ri,R2;o;)Gf*(R2,Ri;a;)] . (5) 

where J2 is the coupling of the second impurity with the electron gas, f{uj) is the Fermi 
function and G^^*(r, v';uj) is the retarded Green function for electrons with spin a. 

In general, to lowest order, the magnetic behavior of a collection of impurities at coordi- 
nates {Ri} is given by the following Hamiltonian: 

^ = E -^^^S^ • (6) 

where the summation is done over all impurity pairs. The exchange parameters J--^ can be 
written as Jij = JiJjC{Ili, Rj) with the correlation function C(Ri, Rj) given by 



c(R..R,)=l/../(.)i,„|Gr(R,.R,;.)Gr(R.R,;.)l 

= -- /" du;Im IC'^Ri, Rj;ujf] (7) 

The second line is obtained after assuming zero temperature, using the relations 
^^^^(Ri, R2;a;) = (Gj''^'^(R2, Ri;a;))* and taking spin independent Green functions, this 
is exact in the absence of spin-orbit coupling and external magnetic fields. Prom here on 
we drop the spin index in the Green functions. The major contribution to the integral of 
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equation (|7j) comes from states near the Fermi energv|2o|. i.e. from Im [G''*^*(Rj, R-,; ei;')^] 
-^o^ sake of avoiding the time consuming integration of Eq. ((Tj) in the numerical 
method, we use the dominant contribution from the Fermi leveh 

C,,(R„ R,) = -^Im [G'-^HR., R,; ^fY] (8) 

In the free surface, the RKKY interaction is oscillatory with a power law decay 
1/ |Rj — Rjl^ . In the presence of an appropriate surrounding corral this interaction can be 



enhanced and focused^. In this work we exploit the fact that the correlation Ce^(Ri,Rj) 
can be controlled by changing the shape of the quantum corral, i.e. the RKKY interaction 
can be controlled by the quantum corral design. 

Also, this correlation function (8) can be shown to be proportional to the perturbation in 
the LDOS at point R2 produced by a non magnetic impurity located at Ri [5P(R2,Ri)]. 
Therefore, to linear order, the search of magnetic mirages in the spin density is the same 
than the search of mirages in P(R2, Ri); extending the scope of this work to the case where 
the perturbation is non magnetic but still localized in space. 



III. THE INVERSE PROBLEM 



Suppose that N impurities can fit in a lattice with spacing a, and we restrict them to a 
square with side /. The total number of possible configurations (including open geometries) 
is the combinatorial n = (//a)^!/((//a)^ — A^)!A^!. For a typical system of = 30, a = 2.5A 
and / = lOOA, n ~ 10^"^. Therefore, it is impossible to study all these configurations within 
reasonable times. Moreover, in the case of quantum corrals, most measurable quantities are 
expected to be non-smooth-multiple- valley functions of the 2N coordinates of the impurities 
due to the appearance and disappearance of resonant states for different configurations 
of the quantum corral (usually close-shaped corrals produce sharp resonances at discrete 
energies Therefore with a simple relaxation scheme, the system may be trapped in local 
minima of the multiple- valley function. Accordingly, more sophisticated non deterministic 
methods are required. 

In spite of this large number of configurations, for most practical applications it is possible 
to avoid covering all the configuration space. At the same time, it is possible to find good 
enough solutions than can be improved systematically. 
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A. Simulated annealing 



The simulated annealing algorithm IJ, |2l| is a quite general method to obtain minima 



of complicated functions of a large number of variables. In order to use a simulated an- 
nealing method we must define first a cost function E{{sLi}). The cost function is simply a 
mathematical expression that quantifies the desired physical properties of the system; being 
minimum when the target is reached and maximum if undesired properties are present. Ac- 
cordingly, the cost functions used in this work depend on one or more correlation functions. 
Several examples will be discussed in the applications below. 

Once the cost function E{{aLi}) is defined and a random initial configuration is chosen, the 



simulated annealing algorithm consists 



in generating random trial steps and accepting 



them with a probability chosen following the Metropolis algorithm: 

f 1 if ^ 

P{AE) = " (9) 

I exp{-AE/T) if AE > 

where AE is the change on the cost function in the trial step and T is a fictitious temper- 
ature. At each trial step the coordinates of the atoms of the corral are allowed to move 
in a triangular lattice representing the Cu(lll) surface. In the annealing algorithm this 
Metropolis dynamics is repeated while T is gradually reducedj^^. This process allows the 
system to sample the space of variables, staying longer times near better global minimum 
as T diminishes. At T = 0, only favorable [AE < 0) steps are allowed to take place but 
at finite T the system has the opportunity to escape from local minima. Eventually, if 
the fictitious temperature decreases slowly enough, the result will be a global or at least 
a good minimum. A particular property of the present system is, however, that for some 
parameters, as the quality of the minimum increases, the amplitude of the barriers grows. 
Therefore, rather than reducing the temperature, we smoothly adjust it continuously so as 
to maintain a target number of accepted changes. This implies that, in our case, as the cost 
function improves, the temperature rises. Finally, instead of reducing the fictitious temper- 
ature, as in the general simulated annealing case, we slowly reduce the target of accepted 
changes jl^. 

The rate of this reduction and how to define a cost function for a particular problem 
is in general a matter of trial and error and is not unique. Several examples for different 
situations will be given. 
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B. Genetic algorithm 



Due to the presence of resonances, the performance of the simulated anneahng decreases 
as better solutions are found or as a more complicated cost functions are introduced. Since 
we are actually looking for a big resonance, we have to consider a complementary method 
to optimize the configurations further. In order to further improve the optimization we 
combined simulated annealing with another Monte Carlo technique which is sometimes more 



robust, a genetic algorithm 



14| . The genetic algorithm is a learning model which derives from 



an analogy with evolution process in nature. 

In our case, for the quantum corral problem, we have designed a genetic algorithm that 
starts with a population of different corral configurations, termed individuals. In order to use 
genetic algorithms, one must define also a cost function (that we choose to be the same one 
used in the simulated annealing) . The proposed cost function is evaluated for each individual 
of the population. The evaluation of each individual determines its probability to survive to 
the next generation. Statistically only the fittest (lowest cost function) are chosen to survive, 
those which are eliminated are replaced by the offspring resulting from mating survivors. 
In our case, we constructed the configuration of the corrals of the offspring as a random 
combination of the configurations of the corrals of the surviving parents. Accordingly, the 
offprings keep geometrical similarities with their parents. In the general case of a genetic 
algorithm, a random mutation Q| is applied to the offspring in order to maintain diversity in 
the population. In our case, we replace this random mutation step by a simulated annealing 
step for each individual described above. We record the configuration of each individual that 
reaches the lowest cost during the simulated annealing step and use it for the next iteration. 
This combined process allows us to move continuously from a pure simulated annealing to 
a pure genetic algorithm depending on the number of individuals in the genetic algorithm 
and temperature and steps in the simulated annealing. Once a new population of quantum 
corrals is generated, the process starts again retaining the fittest and eliminating the rest. 
This combined procedure was repeated until a good solution was found. 

In Table H] we provide the parameters used to obtain the results shown in the Figures 1 
to 6. 
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Number of 


Targeted 


Sa 


GA 


GA 


Cost 


Case 


impurities 


change 


steps 


steps 


population 


function 


1 


60 


1% 


7 X 10^ 


1 


1 


-9.7, 


2 


90 


1% 


5 X 10^ 


1 


1 


-7.0 


3 


120 


1% 


100"^ 


10^ 


50 


-6.2 



'^for each mutation 

TABLE I: Parameters and results for each of the cases studied 
IV. RESULTS 



A. A single mirage case: Alternatives to the ellipse 

Guided by intuition we may think (as many of us did) that the ellipse was the optimum 
shape for a corral to enhance a perturbation of an impurity located at one focus of the 
ellipse into the other focus. This intuitive though originates in an analogy with geomet- 
rical optics 0]. Unfortunately, the analogy is valid only partially because of the following 
reasons: i) for the characteristic length scales of quantum corrals, geometrical optics is not 
completely valid and non trivial effects can arise from multiple scattering and purely quan- 
tum interference ; ii) the impurities can not be freely placed on the surface, there is always 
an underlying lattice (e.g. triangular for Cu(lll)) that restricts the possible locations of 
the corral atoms. This restriction can have non negligible consequences as pointed out in 
Ref. 16[. iii) There is no theoretical proof that the ellipse is the optimal configuration of 
corral atoms, furthermore, the exact shape could depend on the particular phase shift model 
chosen, i.e. on the specific nature of the corral atoms. 

Given these facts, the following question arises: is there any better atomic configuration 
than the elliptical one? 

To test our method in this simple case, the simulated annealing algorithm was instructed 
to maximize the quantity Ce^(Ri,R2), i.e. to project a perturbation from Ri to R2 as 
efficiently as possible. The cost function was chosen to be simply Ei = — Ce^(Ri,R2). 

The reflection symmetry at the x-axis was imposed on the corral impurities to simplify 
the problem. The distance between the spins at Ri and R2 was fixed to 140 A to resemble 
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the experimental conditions. 

Once the optimization procedure is completed, we plot the correlation function of equa- 
tion (7). In Figure n we show the response function and the quantum corral resulting form 
our optimization process [see Table Q] case 1 for parameters]. We see in Fig. Q that the 
minimum found presents some resemblances to the ellipse but new unexpected features also 
appear. First of all, the atoms of the corral tend to accumulate around the location of the 
perturbation and also around the location of the desired mirage. Atoms on the long sides 
have less importance and are placed in a less dense arrangement. This contrasts with the 
evenly spaced atoms of the early experimental setups 0, Q . A second shell of impurities 
appears to be more efficient than a single shell to confine the electron gas and produces 
better mirages. 

For comparison we have studied the family of con-focal ellipses that can be formed with 
the same number of equidistant impurities (being the perturbation at one focus and the 
mirage at the other). We find that the best mirage formed by the family is a factor 3 
smaller than the one obtained with the structure shown in Fig. ^ 

This comparison with the ellipse demonstrates the power of the optimization technique 
we have developed, finding an alternative and better configuration for the formation of single 
quantum mirage. In the rest of this paper we exploit this technique further to study more 
complex problems. 

B. Design of a quantum corral to generate chosen couplings between three spins 

Returning to the problem of constructing physically an arbitrary Hamiltonian [see Eq. 
(6)], suppose now that we want to build a quantum corral such that generates a target 
magnetic Hamiltonian of three impurities: 

-f^targct = — i7Si ■ S2 — jTSi ■ S3 (10) 

i.e. only a selected set of pairs of the three spins are coupled. We arbitrarily choose the 
coordinates of the spins Ri, R2 and R3 to be in the vertexes of an equilateral triangle of 
side d = 121 A. Note that i) Eq. ^ is a particular case of Eq. (6) where J'23 = and 
J'12 = 1/23 = and ii) the target Hamiltonian has a different symmetry than the desired 
positions of the magnetic spins. That is, we require the magnetic couplings on two sides 
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FIG. 1: (Color online), a) Magnetic response (linear order) to a spin located at Ri (x = 70 A, 
y = 0) or quantum mirage formed as a result of an impurity at Ri Distances are reported in 
A. The corral is optimized to enhance the perturbation at point R2 (x = —70 A, y = 0). The 
underlined term in Hamiltonian in the inset shows which magnetic interaction is associated with 
the correlation function displayed, b) Structure of Ce^(Ri,r) below 5% of the largest peak. The 
design of the corral (white and blue dots) has been done with the Monte Carlo algorithm (see Case 
1 in Table HI). 12 




of the triangle but not on the third. In terms of the electronic spin correlation function 
Cep, this system requires Ce^(R2,R3) = and Ce^(Ri,R2) = Ce^(Ri,R3) to be both as 
large as possible in order to enhance the interaction J . We propose a cost function E2 that 
favors these two requirements that is going to be minimized by the Monte Carlo techniques 
described earlier: 

E2 = — Ce^(Rl, R2) — Ce^(Rl, R3) 

+ |Ce^(Ri, R2) — Ce^(Ri, Ra)! 

+|a,(R2,R3)l (11) 

The first two terms in Eq. (jlip ensure that the interactions are enhanced; the third 
term is used to favor J'12 and jTis to be equal; the last term penalizes the interaction 1/23. If 
the positions of the spins Ri, R2 and R3 are selected a priori and fixed, the cost function 
E2 is just a function of the 2N coordinates {a.j} of the atoms forming the quantum corral. 
We also impose here the reflection symmetry along the x axis. The Monte Carlo algorithm 
automatically minimizes this cost function E2 by varying {aj}. The parameters of the 
calculations can be taken from Table |l] case 2. 

The resulting corral is shown in Fig. |21 and El the original requirements are efficiently 
achieved. In Fig. |21 we see that the optimum configuration found for the impurity atoms 
has dense focusing structures near the source of the perturbation and near the mirages. 
Amazingly, some atoms are automatically located in the middle of the corral to split the 
standing waves towards the "targets" at R2 and R3. In Figure 01 we show the effect of a 
magnetic impurity at R3 which produces a large response at Ri but not at R2 [in contrast 
with the case of an impurity located at Ri (compare with Fig. Let us emphasize 

that the corral is the same in Fig. |21 and |S1 the difference is the position of the magnetic 
perturbation. 

Comparison of Figures El and El shows that the RKKY interaction will induce couplings 
as initially designed. The resulting magnetization shows that S2 and S3 are not directly 
coupled while at the same time the interactions between Si and S2 and between Si and S3 
are enhanced. We emphasize that the positions of the spins and the Hamiltonian chosen 
are arbitrary. We believe that we can find quantum corrals for a large number of targeted 
Hamiltonians using the same approach. 
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FIG. 2: (Color online) Quantum corral optimized to generate the spin Hamiltonian of Eq. (|lfl|) . 
Response function Cej^(Ri,r) due to a spin located at Ri (x = 70 A, y = 0). Note the magnetic 
mirages at R2 and R3. The underlined terms in the Hamiltonian in the inset show the magnetic 
interactions associated with the correlation function displayed. Same conventions and symbols as 
in FigH 
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FIG. 3: (Color online) Quantum corral optimized to generate the spin Hamiltonian of Eq. (|lfl|) . 
Magnetic response function Ce^(R3, r) due to a spin located at R3. Same conventions and symbols 
as in Figs, ^and |21 
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C. Four spins Hamiltonians 



In this section we show the resulting design for four spins with a target Hamiltonian: 

-f^targct = — i/Si • S3 — 1/82 • S4 (12) 

where the positions of the spins form a square. 

In order to achieve the target Hamihonian of Eq. p2|) by means of modulating the RKKY 
interaction, we chose the following cost function: 

E3 = — C*e^(Ri, R3) — Ce^(R2, R4) 



Cep(Rl, R3) — Cep(R2, R4)| 
Cep(Rl, R2)| 
CeF(^2, Rs)! 
Cep(R3, 1^4)1 

Ce,(R4,Rl)| (13) 



The first two terms were introduced to increase the interaction between spins at opposite 
vertexes; the role of the third term is to favor identical coupling between the two pairs 
of spins (J7i3 = = <J, note that any difference will increase the cost function); the 
rest of the terms penalize the coupling between consecutive spins (in order to achieve that 
J'12 = J^23 = 1/34 = 1/41 = 0). For this case we impose symmetry along the "x" and "y" 
axes. The method for optimization is similar to the one used in the cases described above 
but in this case the genetic algorithm plays a more important role. The parameters of the 
calculations can be taken from Table H] case 3. 

The results of minimizing Eq. ()13|1 are shown in Fig. |3 and El Figure ID shows the spin 
density generated by a perturbation at Ri and Figure El the spin density when the magnetic 
perturbation is at R2. Note in Fig. jUthat a magnetic mirage appears at R3 {x = —70 A, 
y = 0) but neither at R2 (x = 0, y = 70 A) nor at R4 (x = 0, ?/ = —70 A). Conversely in Fig. 
the magnetic mirage forms at R4 (x = 0, ?/ = —70 A) but neither at Ri (x = 70 A, ?/ = 0) 
nor at R3 {x = 0, y = —70 A). Therefore, the spin Hamiltonian formed by this structure has 
two pairs of spins interacting independently. Each pair of spins do not couple with the other 
pair, even though the interaction is mediated by an electron gas that is shared by both pairs. 
Such couplings have been recently achieved experimentallyj23|], here we present a structure 
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that achieves a similar objective where the mirages are much closer and couplings between 
crossing pairs are specifically avoided. 



D. Optimization for other purposes: optimal corrals to determine phase shifts 

Let us discuss for the moment the limitations of our approach. Our model describes 
the resonances in the densities of states in terms of an scattering theory which depends 
on a single parameter: the phase shift 6. The perturbation in the density of states due 
to single impurity [2^ or an entire corral 19] has been used to measure the characteristic 
phase-shift of a particular type of impurity. We note, however, that such a procedure is not 
very sensitive to the phase shift and that might be the reason of some controversy in the 



literature |19|, l25l |. The origm of this lack of sensitivity comes from the fact that the systems 
fitted were not optimized to be sensitive to the phase shift. For example, dense and closed 
shaped corrals do not produce densities of states that is sensitive to the phase shift of the 
individual atoms of the corral. It is therefore our goal in this section to demonstrate the 
possibility of optimizing a quantum corral for the determination of the phase shift 6. 

We first describe the optimization of a quantum corral geometry with the aim of measur- 
ing 6. Subsequently, we characterize the corral as an experimental device to measure phase 
shifts. 

Let us first concentrate on Figure El that shows a corral optimized to obtain a couple of 
mirages at symmetric points R2 and R3 due to a perturbation at Ri. As we noted in Section 
II Eq. (8), within linear response, magnetic mirages are proportional to finite differences in 
the density of states under a paramagnetic perturbation [i.e. Ce^(R2,Ri) ~ 6V{Tl2,^i)] 
The latter have been measured experimentally. J, |3| 

As the corral in Figure |21 was forced to have reflection symmetry along the x axis, 
5X'(R2,Ri) and 5P(R3,Ri) are identical. However, if we place an additional impurity 
outside the x axis at a point R4, the symmetry of the system is broken. A difference in 
the amplitudes of the mirages at R2 and R3 due to a perturbation at Ri will reflect the 
contribution of this single asymmetric impurity at R4. Intuitively one would guess that 
the difference will be largest if we put R4 at (Ri + R2)/2 in order to interfere with the 
coupling between Ri and R2 but not with the coupling between Ri and R3. However, 
in order to help intuition, we optimized the corral further so as it gives a large difference 
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FIG. 4: Quantum corral optimized to generate the spin Hamiltonian of Eq. (|12j) . Magnetic 
response function Ce^(Ri,r) due to a spin at Ri (x = 70 A, y = 0). Same conventions and 
symbols as in Fig. ^ Compare to Fig. |SJ 
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FIG. 5: (Color online) Quantum corral optimized to generate the spin Hamiltonian of Eq. (|12|) . 
Magnetic response function Cep(R2,r) due to a spin at R2 (x = 0, y = 70 A). Same conventions 
and symbols as in Fig ^ Compare to Fig. 0J 
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= SV(R2, R-i) — SV(R^, Ri) when there is a fix impurity at R4 while the rest of the 
impurities have refiection symmetry along the x axis and we perturbed the system with an 
impurity at Ri. With that goal we design the following cost function: 

E4 = — Ce^(R2, Rl) — Ce^(R3, Rl) + |Ce^(R2, Rs)! 
— |(Ce^(R2, Rl) + Ce^(R2, R4)) 

-(C,,(R3,Ri) + C,,(R3,R4))| xO.5. (14) 

The first part of Eq. (fT^ is the same as Eq. (fTT|) . The second part is maximum when there 
is a big difference in the mirages at Ri and R3 when there is an additional impurity at R4 
(treated in the limit of linear response). The weight factor (0.5) in the last term in Eq. (fT^ 
controls the relative importance of large mirages vs. large differences. 

After the simulated annealing procedure, the corral obtained with the cost function ()14|l 
is similar to the one obtained in Fig. |21 because the cost function and (fT^ are similar 
27[. Let us concentrate on the amplitudes of the mirages (SV) as a function of the phase 
shift 6 of the impurities of the corral. In the inset of Fig. Efa) we show the amplitude of the 
mirage (either at R2 or R3) due to a perturbation at Ri for the symmetric case (without 
the impurity at R4) as a function of the phase shift parameter 6. 5X'(R2,Ri) is measured 
in units of the density of states at the Fermi level of the unperturbed electron gas. We 
see that the amplitude of the mirages grows monotonically with 6 saturating around 2.5 . 
While the magnitude of 6V could be used to determine, by comparison with an experiment, 
the parameter 6, that would be not enough to verify the validity of the model used to 
construct the corral. As for a different model (i.e. a real phase shift) we would find a 
different way to fit the same data: we need at least two measurements to verify the validity 
of a model with a single free parameter. In Figure El (black line) we plot when there 
is an additional impurity at R4 as a function of a common phase shift 6. The value of A^ 
also increases monotonically with 6 but saturates earlier (around 6 = 1). Note that the 
difference between the two mirages can be of the order of 20 % of the density of states 
at the Fermi level. Therefore, we believe that the contribution of a single impurity can 
be measured within experimental resolution, providing alternative data to corroborate the 
present surface scattering model. 

Finally, let us characterize this double mirage structure as an experimental device to 
measure surface scattering cross sections or phase shifts of arbitrary impurities. Let us 
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assume that the surface of Cu is doped mainly with type one atom (lets say Co) but also 
that there are additional minority impurities on the surface available. Let us say that we 
are not even sure what is the phase shift of our majority impurity. The following question 
arises: is this quantum corral made of impurities of unknown phase shift good enough to 
estimate the phase shifts of any other impurity?. The doted lines lines in Figure |B] have been 
drawn assuming that the impurities of the corral have a phase shift 6 of 0.5 1 and 2 as a 
function of the phase shift 6/ of a different the impurity in R4 and a perturbation with phase 
shift (5 at Ri. We see that is quite insensitive to the phase shift of the impurities of the 
corral, suggesting that only the geometry chosen for confinement determines A^. Therefore, 
a single quantum corral device can be used to "measure" the phase shift of every impurity 
in the periodic table in spite of the uncertainly in the knowledge of the phase shift of the 
impurities in the corral. Note finally, that the cost function was build such that the corral 
was very sensitive to small perturbations at R4. Therefore, our proposed "experimental 
device" is very sensitive to small phase shifts and saturates faster for larger ones than the 
amplitude of the mirages (see inset). Alternatively one might have chosen a different cost 
function to selectively measure higher vales of 5f. 



V. DISCUSSION AND CONCLUSIONS 



We have demonstrated the possibility to design theoretically quantum corrals with mul- 
tiple quantum mirages. In this work we have overcome the problem posed by resonances 
and multiple scattering in quantum corrals that give rise to a both a very large configura- 
tion space and a very irregular response function. This response function structure makes 
deterministic algorithms fail to find good minima. The automated design algorithm is use- 
ful to attack the inverse problem specially in cases where the intuition fails or gives poor 
answers. Instead of covering the huge space of all possible configurations, we propose a non 
deterministic method to explore the possible configurations. For the examples covered here, 
10^ steps were enough to optimize the structure. 

The enhancement of magnetic interactions due to electronic confinement together with 
the automatic design of the quantum corrals allows us to tailor structures that generate 
a variety of spin Hamiltonians involving several magnetic neighbors. We demonstrate in 
this work that it is in principle possible to design quantum corrals such as the geometry of 
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FIG. 6: (Color online) = 52?(R2) — SVill^) difference produced by an additional impurity at 
R4 as a function of the phase shift of all impurities (full line), or assuming a fixed value of the 
phase shift of the corral (red 9 = 0.5i, blue 9 = li, magenta 9 = 2i ) but varying the phase shift of 
the impurity at R4. Units in density of state of free electron gas. 

the system produces magnetic images at predetermined points. This opens the possibility 
to actually construct experimentally magnetic Hamiltonians of any sort imposed a priori. 
This could open a path to study experimentally a number of Hamiltonians which have been 
subject of intense theoretical research along the years. In particular one could study one 
dimensional spin Hamiltonians where there are analytical solutions. Such ability to construct 
Hamiltonians might be relevant in the context of quantum computation. We note in passing 
that one of the proposals to build quantum computers is based on the interaction of nuclear 
spins mediated by an electron gasjll. The physical realization of prescribed Hamiltonians is 
a necessary step in quantum computation. The measurement of isolated spins in quantum 
corrals has been experimentally demonstrated. [J, The exchange Heisenberg coupling was 
theoretically shown to be suitable for quantum computation 2]. The combination of these 
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two ingredients and our method to design optimum quantum corrals opens the possibihty 
to actual realizations of quantum computers with spins on a surface. But, in addition, 
Quantum Computation requires the ability to turn interactions on and off P|. In principle 
the interactions between impurities can be turned on and off by changing the position of 
the Fermi level in the surface band with respect to the bulk band. While the control of the 
Fermi level might be difficult in copper surfaces it is much easier in semiconductor quantum 
wells. 
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